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Abstract 

Recent progress of black hole thermodynamics due to Hawking radiation has been applied to 
universal thermodynamics with universe as inhomogeneous LTB model. Hawking-like temperature 
with quantum corrections has been evaluated using HJ method and the semiclassical part of the 
temperature coincides with that formulated in radial null geodesic method. 
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1 Introduction 

Hawking radiation has been extensively discussed since its discovery and different approaches [3J |31 3] have 
been formulated ti study it. Recently there is a lot of attraction to the newly developed semi classical tunnelling 
method [5l[6l[7].The results agree with those of Hamilton- Jacobi (HJ) method [HI, P at least at the semiclassical 
level. In recent past a fully quantum mechanical description has been formulated by Banerjee et al [10 using 
HJ formulation. They have shown logarithmic correction to the entropy functions in leading order to quantum 
formulation. Usually, Hawking radiation is analysed [TTl [T^ [T31 E] on static background space-times where 
there is a global horizon — the event horizon. But dynamical space time, it is not definite to have an event 
horizon. Recently, Hayward et.al. [H] presented a locally defined Hawking temperature for dynamical black 
holes by using tunnelling method. Subsequently, this approach has been used to universal thermodynamics 
particularly to the FRW universe [Ml |T7l [TBI [12] with apparant horizon as the boundary. Also inhomogeneous 
model of the universe in LTB space time has been discussed recently [30] . 

In the present work we study Hawking like radiation from the dynamic horizon in LTB model of the 
universe using both the approaches. In section 2 Hawking like temperature has been evaluated using the 
radial null geodesic method (popularly known as tunnelling method). Then HJ approach has been used to 
formulated quantum prescription of the Hawking like temperature in section 3. An expression for entropy has 
been presented in section 4. At the end there is a brief summary of the present work in the section 5. 



2 Semiclassical Radial Null Geodesic Method : Hawking like Tem- 
perature 

The inhomogeneous spherically symmetric model of the universe is described by the Lemaitre-Tolman-Bondi 
model with metric ansatz 

ds" = -dt" + -^-—dr^ + R^dnl (1) 
1 + /(r) 
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Here {t, r, 9, (f>) is an orthogonal co moving coordinate with 't' the co-moving time corresponding to a co-moving 

A. and — 
dr ■ - df 



observer, R — R{r, t) is the area radius of the spherical surfaces and ' = ^ and . = Jj. The scalar function 



/(r)(> —1) (known as curvature scalar) identifies the space time as 

(i) bounded if -1 < /(r) < 

(ii) marginally bounded, if, /(r) — 
and 

{in) unbounded if /(r) > 0. 

It is possible to decompose the above space time metric into metric ansatz on the surface of the 2-sphere and 
on the 2D hyper surface normal to the 2-sphere as 

ds^ = habcix'^dx^ + R^dnj (2) 

where = d^ag (-1, ^) is the 2D metric normal to the 2-sphere. The trapping horizon, a hyper-surface 
foliated by marginal spheres is defined by |21j 



d+Rr = 0, I.e., Rr = Vl+TW (3) 

where d± = — \/2 =F ^^j^/^''' (?r^ are null vectors normal to the 2-sphere. Thus the trapping horizon 

coincides with apparent horizon as in FRW space time. 

To study the radial null geodesies, we transform our space time coordinates to painleve-type coordinates by 
the transformation r — >■ R{r, t) and the metric ([T]) now becomes 

ds^ = -\l ^ 1 df - 2 dtdR + + R^dnl (4) 

V 1 + /(r) J 1 + f{r) 1 + /(r) ^ V ) 

In this coordinate system the Kodama vector has the form [52] 

= (v/l + /W, 0, 0, O) (5) 
and the associated conserved energy is given by 

aQ 

c. = -v/TT7m^ (6) 

with S, the classical action of a massless particle. So . " can be interpreted as the energy of the particle 
as measured by an observer with the Kodama vector. 

The differential equation for the radial null geodesic (i.e., ds^ = = dft^) can be written as 

^ = 7?±v/rT7w (7) 

where as usual +/— sign stands for outgoing/incoming null geodesic. As for the universal model we consider 
tunnelling from the outside to the inside of the horizon so we consider only the incoming geodesies. According 
to Parikh and Wilczek [5], the tunnelling probability is characterized by the imaginary part of the action 
corresponding to tunnelling of particles through a barrier (i.e., the classically forbidden region) and we have 

ImS = Im piidR = Im dR dp'j^ = Im I dR (8) 

jRirt jRin Jo JRin Jo R 

where to obtain the last equality we have used the Hamiltonian equation namely, 

^ _dH _ dH 
dpR dpii 

In the above equalities we denote by pR the radial momentum of the tunnelling particle, and Rout are 
positions very close to the horizon with the initial position and Rout, a- classical turning point. Now 
substituting the value of 7? from the equation ([7]) into ^ we have 



ImS = Im 



^<"" dR 
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culm 



dR 



Rin 



(9) 



Thus the tunnehing probabihty T ^ exp { — ^/toS*} can be compared with the Bohzmann factor exp { — ^} to 
obtain the temperature of the trapping (apparant) horizon as 



T = 



Im 



dR 



(10) 



This is the semiclassical Hawking like temperature (except for a factor of ^) of the inhomogeneous LTB universe 
for tunnelling of massless particles across the trapping (apparent) horizon . This problem of 'two' factor has 
been pointed out in |10i 1201 and it is mainly of the following two facts : Firstly, Im Jj^ ""^* PR.dR is not canon- 
ically invariant (i.e., not a paper observable) and secondly, the canonically invariant quantity Im <fl^°"* p^dR 

is normally equivalent to 2Im /^°"* PndR, because for ordinary tunnelling the amplitude of tunnelling is same 
for left <^=^ right, but in our case the trapping horizon is an one way membrane and the amplitude for 
left right are not equivalent. However, in the following section, the HJ method is free from this discrep- 
ancy we shall have correct Hawking-like temperature. 



3 Hamiltonian-Jacobi approach : quantum corrections 

Due to non-static nature of the LTB model there is no time like killing vector, but the Kodama vector 



R 



(11) 



which is time like inside the horizon, takes the role of the time like Killing vector as in stationary black hole 
space time. Accordingly, there is a conserved quantity of a particle moving in this space time as 



dSo 
dt 



dt 

where 5*0 is the classical action of a massless particle. So . " 
an observer with the Kodama vector. We strat with the Klein-Gordan(KG) equation 



(12) 



is the energy of the particle as measured by 
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5^(5^''y=55.)^ = 



(13) 



where the scalar field -0 describes a masless scalar particle. As we concentrate our attention to radial geodesies 
only and further due to spherical symmetric nature of the space-time so it is sufficient to consider the above 
KG equation in the 2D hyperplane {t, r) for which the metric hab is given in equation ([2|). Thus the explicit 
form of the above wave equation is 
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l + f{r)\cH 1 d fl + f{r)\ d^p 



2dr 



dr 



Now writing 



i^(r, t) - exp <j -5(r, t) 



we obtain the differential equation for S as 



R" \ d_ l + f{r) \ dl ^ 
1 + f{r) ) dt\ i?'2 ) dt 



(14) 



(15) 



dS_ 

dt 



1 + f{r) 

i?'2 



dS_ 

dr 



+ ih 



dt^ 



l + f{r)\ d'S 



i?'2 



R'- 



dr'^ 2(1 + f{r)) dt 



d fl + fir) \dS ld/ l + f{r) \dS 
i?'2 J 'dt^2dr [ ) 97 



= 



(16) 
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we now try to solve this complicated partial differential equation (p.d.e.) in perturbative approach with 
Planck's constant as the perturbation parameter, unperturbed term as semiclassical approximation and terms 
in different powers of h correspond to different order quantum corrections. So we write 



Sir, t)^So{r, + Sfc?i'-5fc(r, t) 



(17) 

with fc, a positive integer. If we plug this action ansatz into p.d.e. (jl6p and equate equal powers of h on both 
sides then we have the following set of p.d.e s : 



1 + fjr) 



dSo 
dr 



= 



l + /(r)\ d Sod Si 



i?'2 

fir) 



dr dr 



i?'2 



d^Sn 



R'- 



dr^ 2(1 + /(r)) dt 



d fl + fir) \dSo Id 
i?'2 dt 2dr 



l + /(r) \ dSo 
i?'2 / dr 



= 



(18) 



(19) 



dSi 
dt 



and so on. 



^dSodS2 fl + .f{r) 
i?'2 



dt dt 
^d^Si 



dt^ 



dSi 
dr 



l + .fir ) \ d^S 

i?'2 



dr dr j 

d fl + f{r)\ dSi 1 d fl + f{r)\ dSi 



dr^ 2(1 + /(r)) dt 



i?'2 



dt 2dr 



i?'2 



dr 



= 
(20) 



At a glance it seems that different order p.d.e. s are very complicated, but surprisingly there is lot of 
simplifications by using previous equations of the set. Finally, we obtain exactly same p.d.e for 5*0 as well as 
for Sk, k = 0, 1, 2, .... i.e., 



dSkY_ fl + f{r)\ fdSk' ' 
dt 



i?'2 



dr 



= 



(21) 



Thus different order quantum corrections are not independent, they are proportional to the semiclassical action 
5'o. The proportionality constants can be determined from the dimensional analysis as follows : In the choice 
of units G = c = Kb = 1, the plank constant Ti has the dimension Mp [Mp — Planck mass), so Sk has the 
dimension fT^ = M^'^^, where M is identified as the mass of the universe. So we write 



Sk = akM'^'^So. 
where akS are dimensionless constant parameters. 

As a result the series (ITTt now becomes 



(22) 



S{r, t) = So{r, t) 



1 + Sfcttfc 



M2 



(23) 



Thus solution of 5*0 from equation (|T8)) gives the complete solution for the action S. Using (|12p and (ITSt the 
solution for 5*0 can be written in the integral form as 



So{r, t) 



jdt 



dr 



(24) 



where — /+ sign corresponds to ingoing/outgoing scalar particle and the corresponding wave functions are given 
by (using (fT4)) ) 



exp 




ujdt 



R'dr 



{vttj-r} J yrT7{yr+7-^;} 



(25) 
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and 



exp 



wdt 



R'dr 



(26) 



At this point it is to be noted that in crossing the horizon the matrix hab of the metric corresponding to 2D 
hyperplane — (r, t) becomes —hab- i-e., the metric coefficient gu and change their signs, and consequently, 
the above time integration may have imaginary part and hence contribute to the probabihties for the incoming 
and outgoing particles. So the probabilities are given by 



Pin — 



exp 



and 



exp 



2 f h 




Lodt 



ujim 



R'dr 



(27) 



Im 



Lddt 



Lolm 



R'dr 



[VTTJ-r] J VTT7{yTT7-i?} 



(28) 

In the classical domain the horizon is no longer a barrier (absorber in this case) and everything can go out so 
we must have limn^QPout = !> which gives 



Im 



and as a result Pin simplifies to 



exp 



ujdt 



{vttj-r] 



Ldlm 



R'dr 



VTT7{VTTf-R} 



(29) 



R'dr 



VTTf{VTT7-R} 



(30) 



so from the principle of 'detailed balance' we have. 



Pout = exp |-— ^ 

and comparing with ([50)1 we obtain the Hawking like temperature with quantum corrections as 

-1 



(31) 



h ( f h 

T, = - l + S.a. — 



where 



R'dr I I f h 



(32) 



' 4 \ "7 ^{l + f{r)){R-VT+W)} J 

is the semi-classical Hawking-like temperature at the trapping horizon of the LTB model. 

Thus the quantum corrections appear as a multiplicative factor to the semiclassical temperature. Different 
order quantum corrections appear additively as an infinite series. The parameters au are chosen in such a 
manner that the infinite series converges. 



4 Entropy of the horizon with quantum corrections 

We shall now express the entropy of the trapping horizon from the thermodynamical law namely, 

dM = ThdSh 
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where M is identified as the Mishner-Sharp mass [20] 

Now using p2|) in (|33p we have the quantum correction entropy of the LTB universe bounded by the trapping 
horizon as 

Due to inhomogeneity of the space-time, one can not evaluate the integrations, so it is not possible to 
derive the semiclassical Bekenstein-Hawking entropy formula from the above expression. Also we can not 
predict whether the leading order quantum correction gives logarithmic contribution or not. However, the 
above integrand can be simplified to some extend by choosing a'f. = a'^ then the summation term can be written 
in compact form as (l — ^) ^ and hence the form of Sh simplifies to 



^_2/" dR I r R'dr 1 

'"^7 i-f I v/i + /M(^-v/i + /W)J 

It should be mentioned here that such choice of the proportionality parameter aj, was used by Banerjee et 
al [To] in deriving one loop back reaction effects in the space-time for static black holes. 



5 Summary 

In the present work we have formulated a quantum prescription of the Hawking like radiation from the trap- 
ping horizon of the inhomogeneous LTB model of the universe. Due to non-static nature of the space time 
we have used Kodama observer and the energy of the particle is associated with the corresponding invariant 
energy. The outgoing and incoming probabilities are matched using the principle of detailed balance. The 
semiclassical part of the quantum corrected temperature using HJ approach is compared with that obtained 
from tunnelling approach where there is a discrepancy of a factor " two" . Finally, an expression for the horizon 
entropy has been presented and it is not possible to deduce the standard entropy-area formula at the semi clas- 
sical level due to its complicated form. For future work, we shall attempt to resolve the factor two discrepancy 
for non-static metric and try to find any physical interpretation to the different order quantum correction terms. 

Acknowledgement : 

RB and NM want to thank West Bengal State Govt, and CSIR, India respectively for awarding JRF. Authors 
are thankful to lUCAA, Pune as this work was done there during a visit. 

References 

[1] Hawking, S. W -.-Nature 30 248(1974); Commun. Math. Phys. 43, 199(1975). 

[2] Hartle, J.B., Hawking, S. W. :- Phys. Rev. D 13, 2188(1976). 

[3] Gibbons, G.W., Hawking, S.W. -.-Phys.Rev.D 15, 2752(1977). 

[4] Christensen, S. M., Fulling, S. A. :- Phys. Rev. D 15, 2088(1977). 

[5] Parikh, M. K., Wilczek, F. :- Phys. Rev. Lett. 85, 5042(2000). 

[6] Parikh, M. K. :- preprint 0402166[hep-th]. 

[7] Medved, A.J.M. :- Phys. Rev. D 66, 124009(2002). 

[8] Srinivasan, K., Padmanabhan, T. :- Phys. Rev. D 60, 024007(1999). 

[9] Shankaranarayanan, S., Srinivasan, K., Padmanabhan, T. :- Class. Quantum. Grav. 19, 2671(2002). 



6 



[10] Banerjee, R., Majhi, B.R. :- JHEP 06, 095 (2008). 
[11] Hu, Y.-R, Zhang, J.-Y., Zhao, Z. :- arXiv 0601018. 
[12] Sarkar, S. , Kothawala, D. :- Rhys. Lett. B. 659, 683(2008). 
[13] Rilling, T. :- Phys. Lett. B 660 402(2008). 

[14] Zhang, B., Cai, Q. Y., Zhan, M. S. :Phys. Lett. B. 665, 260(2008). 

[15] Hayward, S. A., Dicriscienzo, Di. R., Vanzo, L., Nadalini, M., Zerbini, S. :- Class. Quantum. Grav., 26 
062001,(2009). 

[16] Cai, R. G., Cao, Li-Ming., Hu, Ya-Peng. :- Class. Quant. Crav. 26, 155018(2009). 
[17] Zhu, T., Ren, JI, R., Singleton, D. :- IJMPD 19, 159(2010). 

[18] Jiang, Ke-Xia., Ke, San-Min, Peng, Dan-Tao, Feng, J. :-preprint hep-th 0812.3006. 
[19] Chowdhury, B. D. :- Pramana 70, 593(2008). 

[20] Biswas, R., Mazunider,N., Chakraborty.S. :- preprint 1106.2286 [gr-qcj. 
[21] Hayward, S. A. :- Phys. Rev. D 53, 1938(1996). 
[22] Kodama, H. :- Prog. Theor. Phys. 63, 1217(1980). 



7 



